Nonlinear localized modes at phase-slip defects in waveguide arrays 
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We study light localization at a phase-slip defect created by two semi-infinite 
mismatched identical arrays of coupled optical waveguides. We demonstrate 
that the nonlinear defect modes possess the specific properties of both nonlinear 
surface modes and discrete solitons. We analyze stability of the localized modes 
and their generation in both linear and nonlinear regimes. 
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The study of nonlinear dynamics in discrete systems 
has attracted a special attention recently due to novel 
physics and possible interesting applications^. In partic- 
ular, it is well known that discrete photonic systems can 
support different types of spatially localized states in the 
form of discrete solitons^^^. These solitons can be con- 
trolled by the insertion of suitable defects in an array, 
as was suggested theoretically^'^ and also verified exper- 
imentally for arrays of optical waveguides^. Defects may 
provide an additional physical mechanism for light con- 
finement, and they can support both linear and nonlinear 
localized modes, which has been studied theoretically for 
different nonlinear models^^^^^ and observed experimen- 
tally in one-dimensional photonic lattices^. 

In this Letter, we introduce a novel type of nonlin- 
ear defect in waveguide arrays closely linked with the 
recently discussed phase-slip defects in two-dimensional 
photonic crystal^iii. In particular, we demonstrate 
that two semi-infinite mismatched identical arrays of op- 
tical waveguides can support a variety of linear and non- 
linear localized modes with the specific properties of both 
discrete soliton^ and nonlinear surface modesi^^i^. We 
analyze stability of the localized modes and their gener- 
ation in both linear and nonlinear regimes. 

We consider an array of nonlinear optical waveguides 
created by two semi-infinite identical mismatched arrays, 
as shown in Fig. 1. In this array, two mismatched waveg- 
uides at the sites n = m and n = m-\-l interact with a dif- 
ferent coupling parameter, V' so that the coupled- 
mode system can be described by the discrete equations 
for the normalized mode amplitudes En^ 
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where En are defined in terms of the actual electric fields 
En as En = {2Xor]o/7rnon2)^^'^En, where Aq is the free- 
space wavelength, 770 is the free-space impedance, 712 



and no are nonlinear and linear refractive indices of each 
waveguide, and 7 (±1) defines the type of nonlinear ity. 
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Fig. 1. Transverse profile of the refractive index for an ar- 
ray of weakly coupled nonlinear optical waveguides with 
a phase-slip defect located between the waveguides at the 
sites m = and m = 1. Top: > V ^ bottom: V < V. 

We look for stationary solutions of Eqs. ([l]), ([2]) in the 
form En{z) = Ene-x.-p{i/3z)^ and consider first the linear 
waveguide array (or the limit of low beam powers) when 
7 = 0. In this case, we expect that localized modes may 
exist for V > V only, since decreasing the ratio V jV 
decouples the chain into two pieces, and each of the semi- 
infinite chains does not support surface mode^ii^. We 
search for localized solutions of the form En = 
for n < and En = B^^'^~'^~^\ for n > m + 1, where 
1^1 < 1. After some algebra, we obtain /3 = V[£, -\- (1/0] 7 
B/A = ^{V'/V) and ^ = ±\V/V'\ so that indeed, local- 
ized modes require the condition V > V. 

For ^ = |V^/y| and A = B we obtain unstaggered 
localized modes [see Fig. 2(a)]: En = A|l//V'|l^-^l , for 
n < m, and En = A\V/V'\''-'^-^ , for n > m + 1 with the 
propagation constant P/V = \V/V'\ -\- \V' /V\. Similarly, 
for ^ = —V/V^ and A = —B we obtain the corresponding 
staggered localized modes [see Fig. 2(b)]. 

Next, we consider the nonlinear case described by 
the stationary form of Eqs. ([l]), ([2]) for 7 > 0. For a 
given value of /3, the system of the stationary equations 
is solved numerically by a multidimensional Newton- 
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Fig. 2. Linear unstaggered (a) and staggered (b) localized 
modes at the phase-slip defect for V /V = 1.3 {N = 
100, m = 50). 



Raphson scheme. As we are interested in the modes lo- 
calized near the defect, we look for the modes with the 
mode maxima near the slip boundary that decay quickly 
away from the bond impurity. We find that the results 
vary depending on whether V > V or V < V. 
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Fig. 3. Power vs. propagation constant for several fami- 
lies of nonlinear localized modes at the phase-slip defect 
for VyV = 1.3. Solid (dashed) curves denote stable (un- 
stable) branches. Inset: Minimum power to destabilize 
the fundamental mode vs. coupling mismatch. 

Figures [3] and [4] depict the mode power vs. propagation 
constant and show specific examples of the mode profiles, 
for the case V > V. First, one of the nonlinear modes 
extends all the way down to zero power, and it generalizes 
the linear mode found above [see the curve (a-c) in Fig. [3] 
and Figs. [4l(a,c)]. This mode becomes unstable above a 
certain threshold power, and it transforms into an odd 
mode centered at any of the two equivalent sites coupled 
by the bond impurity [see Fig. [^b)]. This result can be 



easily understood from the known instability of an even 
mode for a discrete homogeneous lattice: As the power 
is increased, the effective coupling is decreased and the 
distinction between V and V becomes blurred. In the 
high-power limit, the even-like localized state "feels" as 
inside an homogeneous lattice; hence the onset of insta- 
bility. The inset in Fig. [3] shows the critical power needed 
to destabilize the fundamental, even-like localized mode. 
The most interesting feature of this curve is that it rises 
very steeply as soon as the ratio V /V is slightly above 
one, followed by a slow, almost linear-like growth. This 
suggests that a tiny amount of mismatch is enough to 
stabilize the even-like mode at the phase slip. 
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Fig. 4. (a-e) Examples of the nonlinear localized modes 
marked by the letters a,b,c,d, and e in Fig. [3l 

For larger powers, we find novel types of nonlinear 
modes localized at the phase-slip defect [see Figs. |4l^d,e)] 
resembling a bound state of two simpler modes. These 
modes resemble the so-called twisted modes found earlier 
in the homogeneous chain^^^^, and they exist only above 
a certain power threshold. The complementary unstable 
mode looks like the twisted mode with "shoulders" [See 
Fig. 4(e)]. 

In the case of a weaker bond defect, i.e. V < V ^ lin- 
ear localized modes do not exist. This result is consistent 
with the case of surface mode^*^ where a certain power 
threshold is required to support a mode localized at the 
edge of the waveguide array. Similarly, in this model the 
localized mode appear for finite powers (see the branch 
marked 'a' in Fig. [5|). The mode profiles are similar to 
those presented in Figs. [H^ a-e). In this case, all localized 
modes are strictly nonlinear i.e., they disappear in the 
limit 7^0. As a result, most of those modes are un- 
stable, and only two modes corresponding to the lower 
power are stable. 

We should mention that all modes for 7 > remain 
unstaggered, and staggered modes in this model appear 
only for 7 = —1, being found through a simple trans- 
formation En (— l)^£^n applied to all types of modes 
discussed above. 



Fig. 5. Power vs. propagation constant for several local- 
ized states at the phase-slip defect for V jV = 0.7. Solid 
(dashed) lines denote stable (unstable) branches. 

Finally, we analyze the generation of the defect modes 
by an input beam sent to one of the waveguides of the 
phase-slip defect. For weaker coupling {V < V) and 
low powers, we observe no power localization near the 
defect sites, and the input power diffracts as in homoge- 
neous arrays [see Fig. [6l^a)]. However, when we increase 
the input power, we are able to generate the asymmet- 
ric nonlinear defect mode that correspond to the lowest 
branch ( marked with 'b') in Fig. [5] and finite powers, 
similar to the excitation of discrete surface soliton s^^i^^ . 
On the contrary, the surface mode is always generated for 
the case of stronger coupling, V > V ^ when the defect 
mode exists in the linear regime, as shown in Figs.[6l[c,d). 

In conclusion, we have introduced and described novel 
types of nonlinear defect modes localized at a phase- 
slip defect in an array of nonlinear optical waveguides. 
We have demonstrated that these localized modes pos- 
sess many specific properties of both discrete solitons and 
nonlinear surface modes, and they can be easily gener- 
ated in both linear and nonlinear regimes. 
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Fig. 6. One-waveguide excitation of the defect modes. 
Top: V /V = 0.7, i.e. for a weak bond defect, for the in- 
put amplitude (a) Vq = 1.0 and (b) Vq = 2.0 (b). Bottom: 
same as above but for stronger bond defect, /V = 1.3. 
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